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ABSTRACT 



We present slitless spectra of the Narrow Line Region (NLR) in NGC4151 
from the Space Telescope Imaging Spectrograph (STIS) on HST, and investigate 
the kinematics and physical conditions of the emission line clouds in this 
region. Using medium resolution (~0.5A) slitless spectra at two roll angles 
and narrow band undispersed images, we have mapped the NLR velocity field 
from L2kpc to within 13 pc (Ho = 75kms^^Mpc~^) of the nucleus. The inner 
biconical cloud distribution exhibits recessional velocities relative to the nucleus 
to the NE and approaching velocities to the SW of the nucleus. We find 
evidence for at least two kinematic components in the NLR. One kinematic 
component is characterized by Low Velocities and Low Velocity Dispersions 
(LVLVD clouds: |v| <400kms^^, and Av< 130kms^^). This population 
extends through the NLR and their observed kinematics may be gravitationally 
associated with the host galaxy. Another component is characterized by High 
Velocities and High Velocity Dispersions (HVHVD clouds: 400 < |v| ^1700 
kms^^, Av> ISOkms^^). This set of clouds is located within LI" (~70pc) of 
the nucleus and has radial velocities which are too high to be gravitational in 
origin, but show no strong correlation between velocity or velocity dispersion 
and the position of the radio knots. Outflow scenarios will be discussed as the 
driving mechanism for these HVHVD clouds. We also find clouds characterized 
by Low Velocities and High Velocity Dispersions (LVHVD clouds: |v| <400 
kms"^, Av> ISOkms"^). These clouds are located within 3.2" (~200pc) of 
the nucleus. It is not clear if the LVHVD clouds are HVHVD clouds whose low 
velocities are the results of projection effects. 

Within 3.2" (~200pc) of the nucleus, the [0III]/H/9 ratio declines roughly 
linearly for both the High Velocity Dispersion (HVD) and LVLVD clouds. Since 
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the ionization parameter is proportional to r^^n^^, it appears that the density, 
n, must decrease as ~r^^ for the clouds within the inner ~3.2". At distances 
further from the nucleus, the [0III]/H/3 ratio is roughly constant. 
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1. Introduction 

Current models describing Seyfert galaxies consist of an active nucleus surrounded by 
a dense (~ 10^ — lO^^cm^^) broad emission line region (BLR) which lies well within ~ 1 pc. 
Extending to larger radii (~ 1 — 1000 pc) is a lower density (^ lO^cm"'^) region where the 
narrow permitted and forbidden lines of the NLR are generated. Nuclear activity is believed 
to be generated by mass infall from a hot accretion disk surrounding a massive black hole. 
X-rays released by this process irradiate the accretion disk generating a biconical wind. 
This nuclear wind may then be responsible for the outward acceleration of clouds into the 
NLR. 

Radio emitting jets have also been observed in Seyfert galaxies. The jet is believed to 
be aligned with the rotation axis of the accretion disk and may provide another mechanism 
for cloud acceleration. Line profile studies of Seyferts with multi-knot radio jets have 
shown that substructure in the [OIII] emission line profiles is correlated with the spatial 
locations of radio jet knots, indicating that the interaction with the radio plasma may 
produce density enhancements in the NLR gas ( [Whittle et al. 1988| , |Aoki et al. 1998|) . By 
sweeping up gas the radio plasma may also serve to provide a path for nuclear photons, 
thus facilitating ionization in the extended emission line region. 

At distances large compared to the accretion disk the gravitational influence of the host 
galaxy dominates the velocity field. Gas in this region is characterized by orbital rather 
than radial or turbulent velocities. 

Although this complex model has emerged, it is far from verified in every detail and 
the interaction between the NLR and its local environment is still not well understood. 
Emission line profiles show abrupt spatial variations, thus supporting the theory that 
multiple processes are influencing the NLR gas. The blue asymmetry of the narrow line 
profiles has led various authors to advocate either infall ( [DeRobertis and Shaw 1990| ) or 
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outflow (e.g. [Schulz 1990| , [Whittle, 1992| ). Line profiles have been modelled and matched 
to observations for wind models which confine and drive the narrow line clouds ( [Krolik 
et al. 19811 , [Krolik and Vrtilek 1984| , [Schiano 1986aj, [Schiano 1986b[ ). Correlations have 
been found between substructure in the emission line profiles and the position of the radio 
jets ( [Whittle et al. 19"88[ , [Winge et al. 19991) . However, distinguishing among the various 
dynamic processes is inhibited by the low spatial resolution or limited spatial mapping of 
the NLR. Complete, high spatial resolution mapping of the complex NLR gas distribution 
will be an important contribution in distinguishing among the models. 

Ground-based studies of the NLR have been limited to spatial resolutions of ^ 50 pc 
for the closest Seyfert galaxies. In general they have relied upon long-slit spatial line profile 
studies to understand the dynamics of the NLR and its dependence upon the BLR ( [Schulz] 



1990[ , [Veilleux 1991[) . Fabry-Perot observations have mapped the entire spatially extended 



NLR, but these observations are also ultimately limited to coarse spatial resolution by 
seeing conditions (~60 pc resolution for NGC 1068, [Cecil et al. 19901 ). 



With the advent of HST, high spatial resolution imaging of Seyfert galaxies in the 
optical has exposed a rich, complex and varied NLR. Narrow band WFPC images of [OIII] 
and Ha + [Nil] have revealed the morphological structure of the NLR for a large sample 



of Seyfert galaxies (NGC 4151: Evans et al. 1993 , Boksenberg et al. 1995 ). Long slit 
observations have provided kinematic information at a few select positions within the NLR 
(e.g., NGC 4151: [Winge et al. TW7[ , [Winge et al. im% Mrk3: [Capetti et bTW^ ). While 
this complement of observations has expanded our knowledge, the lack of full kinematic 
mapping still inhibits our understanding of the environmental forces infiuencing the NLR. 



As one of the brightest and closest {yrecessionai ~997kms~^, [Pedlar et al. 199^ ) 
Seyfert 1 galaxies, NGC 4151 is an ideal candidate for studying the full kinematic field of 
the NLR emission line gas. Now, with STIS slitless spectroscopy the entire dispersed NLR 
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can be imaged with a spectral resolution of 0.276 A pix^^ at [OIII] and a spatial resolution 
of 0.051" pix-i. 

Our previous STIS slitless spectroscopy presented a kinematic map of the NGC4151 
NLR (|Hutchings et al. 1998| ). In this paper we present a more complete kinematic 



map of the NGC 4151 NLR at a second roll angle with expanded wavelength coverage 
encompassing three wavelength ranges (~ 3700 - 3800 A, ~ 4800 - 5100 A, ~ 6300 - 6850 A) 
to include emission lines with a range of ionization potentials. Expanding upon our 
earlier measurements (|Hutchings et al. 19981) we present a table of velocities and velocity 



dispersions for 57 emission line clouds. [OIII] and H/5 fluxes have been measured for 34 
of the clouds in this sample. Where possible [Oil], [01], and [SII] fluxes have also been 
measured. The velocities and fluxes have been employed to search for radial structure in the 
kinematic and ionization distributions. Correlations with the spatial location of the radio 
jet have been investigated. Outflow, infall, and gravitational models are also examined 
within the context of the velocity, velocity dispersion, and flux distributions. 

The observations are discussed in Section |], the image processing is described in 
Section ^, the velocity measurements in Section ^, and the flux measurements in Section ^. 
The results are presented in Section 6 and discussed in Section 7. 

2. Observations 

Slitless spectra of the nuclear region of NGC 4151 were obtained using the Space 
Telescope Imaging Spectrograph ([Kimble et al. 1998| , [Woodgate et al. 1998|) aboard the 



Hubble Space Telescope. These CCD observations were acquired at a single epoch, with 
the dispersion direction oriented -45° E of N. Three medium dispersion grating tilts were 
selected to permit observation of the spectral region encompassing the [Oil] A3727 emission 
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(Figure |l]), the H/3 and [OIII] A5007 emission (Figure |) and the [01] A6300, Ha+ [Nil], 
and [SII] emission structure (Figure ^. Each format spans a 50" field-of-view in the 
spatial direction. The following emission lines were also detected: HI 2 A3751, [FeVII] 
A3760, Hll A3771, HIO A3798, H9 A3835, [Nelll] A3869, H8 + Hel A3889, and He A3970 
+ [Nelll] A3967. For these observations the resolution is defined by the spatial extent of 
the emission feature. Assuming a 2.0 pixel emission feature, which is equivalent to the 
average instrument resolution for a point source, the resolution is 0.55 A (45 km s~^) at 
[OH] A3727, 0.55 A (33 km s~^) at [OIII] A5007, and ~l.lA (51 km s"^) at Ha. However, 
the near-nuclear emission line clouds typically have velocity broadened widths which are 
then convolved with the Optical Telescope Assembly (OTA) plus instrument point spread 
function (PSF) to produce observed widths of J^9 pixels FWHM, rendering it difficult to 
deconvolve the Ha + [Nil] emission line clouds. 

To assist in identification of the dispersed cloud structure and to facilitate measurement 
of the cloud velocities and velocity dispersion, observations were also acquired with the 
STIS narrow band [OH] fiher and the WFPC2 narrow band filters: F502N, F658N, F390N, 
F673N, and F656N with the Wide Field (WF) camera. The F502N image was used as the 
primary undispersed image. Its FWHM bandpass encompasses emission line clouds with 
velocities extending from — 1700 km s~^ through ~500kms"^ with respect to the systemic 
velocity. The broadest wavelength coverage was provided by the STIS [Oil] filter whose 
bandpass includes velocities extending from — 2700 km s~^ through 2350 km s~^ with respect 
to the nucleus. The STIS [Oil] filter was used to confirm cloud identifications and the STIS 
dispersed images at PA = 104° and PA = 45° were used to search for high velocity clouds not 
contained in the F502N [OIII] bandpass. With the exception of the F502N [OIII] image, the 
WF narrow-band filter observations were acquired at a different epoch than the dispersed 
images. Consequently, the position angle for the WF images is at -41° E of N, whereas the 
position angle is 45° E of N for the dispersed images and the WF [OIII] image. 



Previous STIS data ([Hutchings et al. 19981) were acquired using the same instrument 



configuration as the current [OIII] dispersed data. The data acquired at this earher epoch 
were oriented with the dispersion direction 14° E of N. A log of the observations is presented 
in Table 1; entries marked with an " represent the earlier data. The archival WFPC2 data 
were acquired with the Planetary Camera (PC) channel. The higher resolution of the 
PC compared to the WF camera (0.0455 "/pixel versus 0.996 "/pixel) is advantageous for 
identifying the numerous emission line clouds in the near-nuclear region. These clouds are 
annotated in Figures |^, |^, and ^. 

3. Image processing 

Minimal post-pipeline processing was required for the WF and PC images used to 
determine the spatial location of the [OIII] emission line clouds relative to the NGC 4151 
nucleus. This additional processing consisted of rescaling the images to the STIS spatial 
scale (0.05071 "/pixel) and then rotating the WF and PC images to the same orientation 
as the STIS image. 

For the STIS slitless spectra, scattered continuum light and the host galaxy background 
were subtracted from the post-pipeline processed images. This was done to enhance the 
contrast of low surface brightness clouds for the velocity measurements and to reduce the 
background for the flux measurements. Our background subtracted image exhibits no 
evidence for uncorrected absorption of Ha near the nucleus. 

The background template was constructed by median flltering the image along the 
dispersion direction. This template flt the background well at distances greater than 1.6" 
from the nucleus. However, the inner wings and the core of this spatial template varied 
along the dispersion direction, due at least in part, to the variation of the continuum as 
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a function of wavelength and the scattered hght from the nucleus and near-nuclear clouds 
in the emission lines. Consequently, multiple templates were constructed as a function of 
wavelength to subtract the near-nuclear region. Far-field template construction consists 
of masking the emission line regions, and median filtering the image along the dispersion 
direction. The near-nuclear region (±1.6") is then masked and the far-field template 
subtracted from the entire image. The near-nuclear templates are similarly constructed, 
using rows within ±1.6" of the nucleus, but typically only 30 columns are averaged to 
preserve the spectral signature of the continuum. For the emission line regions, the 
near-nuclear background template is an average of the neighboring templates redward and 
blueward of each emission line region. 

4. Velocity Measurements 

The [OIII] A5007 emission is bright and unblended with other emission lines, and 
therefore the [OIII] dispersed image and the F502N PC images were selected as the 
primary images for determining the cloud positions and dispersions. Since the NLR 
presents a well-ordered velocity field, associating the dispersed clouds with the undispersed 
counterparts is not difficult in general. In congested regions where these identifications are 
more difficult, measurements of the spatial locations of the clouds in the WFPC2 and STIS 
images by various co-authors were required to agree within 1 pixel. 

Multiple approaches were employed to measure the velocity displacement and dispersion 
of the emission line clouds. Measurements were made using both background subtracted 
and non-background subtracted images, using automated and purely interactive cloud 
identifications and fits. Checks were also performed using the STIS [Oil] undispersed image 
with the H/3 and [01] dispersed images. The errors quoted in Table 2 are the one sigma 
spread in the measurements obtained by multiple co-authors using multiple techniques. 
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The quoted errors dominate the statistical errors associated with Gaussian fits to the fine 
profiles. The narrow band filter data were inspected to check for ghosts and the presence of 
the clouds in the various emission lines. 

All cloud velocities were determined relative to the NGC4151 nucleus. The location of 
the nucleus was determined by fitting a Gaussian to the continuum of the non-background 
subtracted image in both the spatial and dispersion directions. Positions and velocity 
dispersions for the individual clouds were determined by fitting Gaussians to the extended 
spectral region encompassing the clouds. In general, the extent and proximity of the clouds 
to one another required fitting multiple Gaussians with a linear background to an extended 
region neighboring the cloud. With this method, the cloud centers with respect to the 
nucleus and the cloud FWHM are determined for both the dispersed and undispersed 
images. Both the relative position and FWHM of each cloud were then corrected for 
anamorphic magnification. This correction demagnifies the STIS pixel size in the dispersion 
direction by 11%. The velocity dispersion is then determined by subtracting the instrument 
profile convolved with the spatial extent of the cloud, defined as the FWHM of the 
undispersed cloud, in quadrature from the dispersed FWHM of the cloud. 

At large distances perpendicular to the dispersion direction, a small misalignment in 
the rotation of the STIS and PC images translates into an error in the cloud velocity. At 
the location of cloud 44, which lies 10.9" from the nucleus, a 0.3° error in the alignment of 
the STIS and PC images corresponds to a velocity shift of 20kms"^ in the [OIII] emission. 
Consequently, the outer NLR cloud velocities are small and potentially consistent with zero. 

In addition to the manual method, an automated procedure was written for matching 
the dispersed emission line features with their counterparts in the narrow-band images and 
calculating the radial velocities and velocity dispersions for the individual clouds. With this 
method, the images are aligned and registered. Then "slices" parallel to the dispersion axis. 
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at the same spatial offset from the zero-point, are plotted for each image. Using a graphics 
cursor, features are marked and labeled interactively such that peaks in the dispersed 
image slice are matched with peaks in the undispersed image slice. Widths and centers of 
profile cuts at 80%, 50%, and 20% of the peak height are evaluated, similar to the standard 
techniques used in analysis of emission line profiles (e.g. [Whittle 1985|) . The relative shift 



between the centers obtained in any of these ways is converted into a shift in wavelength 
and then a velocity. The velocity dispersion is determined by subtracting the width of the 
undispersed image in quadrature from the width in the dispersed image, converting to a 
width in angstroms and then to a velocity width. 

The aforementioned methods were suitable for associating the majority of the dispersed 
clouds with their undispersed counterparts. However, given the cloud congestion in selected 
near-nuclear regions and the potential for spatial-spectral ambiguity in the dispersion 
direction with slitless spectra, a third method was employed. This method was used to 
confirm cloud identification/positions in the dispersed image (i.e. cloud regions 7-10-12 
and 19-22-23), as well as identifying cloud velocities and positions for those clouds which 
were either too faint to be seen in the PC image or redshifted out of the PC F502N filter 
bandpass (i.e. clouds with velocities ^500 kms~^). By employing dispersed images at two 
different position angles the velocity and rest position of clouds with respect to the nucleus 
can be determined. Figure |^ illustrates the method. Here we see that the second image has 
been rotated to the position angle of the first image and the images have been aligned at 
the nucleus. The same cloud is identified in each image (i.e. cloud Al, A2). To determine 
the undispersed cloud position, a line is extended from the cloud parallel to the continuum 
in each image. The intersection of these two lines is the location of the undispersed cloud 
(i.e. point A). The distance from the undispersed to the dispersed position should be the 
same in both images, i.e. A-Al and A-A2 should be equal. Requiring that the velocity shift 
be the same in both images is a powerful tool for resolving cloud identification ambiguities 
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in congested regions. 

Since several of these clouds are faint, a calculation was used to verify the cloud 
identifications determined from the geometric method described above. For the analytic 
method, the position {xmr,Jmr) of the cloud with respect to the nucleus in the archival 
(unrotated) image was measured, and then the ansatz y-position (y^) of the cloud in the 
rotated image was inserted into the combined equations 1 and 2 to determine the x-position 
of the cloud in the rotated image, (x^ + A). Then the velocity and spatial position of the 
cloud in the rotated cycle 7 image were calculated. 



(x; + A) = Xmr + A(l - cos{9)) 



y'o = Vmr + ^sin{e) (2) 



Xmr = {xo + A)cos(^) + yoSin{9) (3) 



Vmr = -{xo + A)sin{9) + y„cos{d) (4) 

In the rotated cycle 7 frame the spatial coordinates of the cloud with respect to the 
nucleus are (x'^, y^), in the unrotated archival frame the coordinates with respect to the 
nucleus are (xo,yo), is the angle between the position angles of the two images, and A 
is the x-distance from the spatial location to the dispersed location of the cloud. Clouds 
whose velocities were confirmed analytically are presented in Table 3. 

A velocity dispersion was calculated for each cloud. The instrument-convolved spatial 
width (A Vi^pc) was determined by a Gaussian fit to the cloud in the F502N PC image. 
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This width was subtracted in quadrature from the Gaussian width fit to the cloud in the 



dispersed STIS image (Avi^sris)- The velocity dispersions {Avi = JAvf^gj,jg — Avf^pQ) 
are presented in Table 2. Systematic errors dominate the statistical errors. The errors 
quoted in the table are the 1 a spread of the measurements by the various authors. 

For the high velocity clouds presented in Table 3, there are no corresponding clouds 
identified in the PC image because these are relatively faint clouds which have been Doppler 
shifted to to the edge or outside the F502N filter bandpass. For these clouds, the velocity 
dispersion was calculated by subtracting in quadrature the average instrument convolved 
spatial width (< Avpc >«) for a set of clouds who have approximately the same value for 
A Vi^sTis ^^ the high velocity cloud. 

In order to determine whether the high velocity, high velocity dispersion clouds might 
simply be the superposition of several low velocity dispersion clouds, we performed a Monte 
Carlo simulation. An ensemble of 3 or 4 clouds with a Gaussian FWHM of a typical LVD 
cloud were superposed with fluxes ranging over a factor of 4 and velocities separations of 
up to 300kms~^. The observed distribution of clouds on the sky makes this a reasonable 
number to superpose, since by symmetry a cylinder of typical cloud diameter, oriented 
perpendicular to the bicone axis, encompasses 4 clouds. In addition, fluxes fainter than 
about 1/4 of the brightest in a group have no visible effect on the combined profile. In over 
95% of the cases, the resultant blends are very structured or multiple-peaked. The observed 
high velocity dispersion profiles are all smooth. Thus, we conclude that the broader line 
clouds do indeed have intrinsic high velocity dispersions. 
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5. Flux Measurements 

H/? and [OIII] fluxes have been measured for 34 of the 57 tabulated clouds. In addition, 
we have measured the [Oil], [01], and [SII] flux for a subset of bright, relatively isolated 
clouds. 

Flux measurements were executed using the background subtracted images. Initially, 
clouds were flt with a linear background and the minimum number of Gaussians, in the 
dispersion direction, required to obtain a reasonable fit to the cloud and its neighbors. In 
general, neighboring clouds were fit simultaneously with the target cloud because the wings 
of adjacent clouds overlap. Fits were performed in sequential rows spanning the fiux peak in 
the spatial direction until the cloud was indistinguishable from the background. The total 
flux for each cloud was determined by summing the area associated with the Gaussian flt to 
each row encompassed by the cloud. Typically clouds spanned ~5— 7 rows (~0.25 - 0.35") 
in the spatial direction. Flux ratios are calculated using the same number of rows for the 
flux of each species. The fluxes listed in Table 5 for H/9 and [OIII] may have lower values 
than the fluxes listed in Table 4 because fewer rows were summed to obtain the [Oil], [01], 
and [SII] fluxes, and hence the [OIII] and H/3 fluxes and flux ratios presented in Table 5. 

Gaussian flts to the clouds tend to underestimate the wings of the line proflles. A 
sub-sample of relatively isolated clouds with sufficient signal-to-noise in the [Oil] image 
were modelled using both Lorentzian and Gaussian line profiles. In general, the Lorentzian 
provides a better fit to the wings. However, the line profile wings are probably elevated 
by scattered light within the emission line region. Due to the complexity of the emission 
structure it is difficult to determine this background. Consequently, the difference between 
fiuxes obtained with the Gaussian and Lorentzian line profile models is treated as a measure 
of our systematic error. The errors quoted in Table 5 are the difference between the fluxes 
calculated from the Lorentzian and Gaussian flts to the data. These errors are positive 
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indicating that the Lorentzian fits yield higher fiuxes, as expected. 

Fits to the [Oil] doublet AA3729,3726 were impeded by the lower fiux of the [Oil] 
emission and the cloud congestion in the narrow line region (NLR) which was exacerbated 
by the 10 pixel (2.8 A) separation of the doublet. Several strategies were employed to 
maximize the number of [Oil] clouds fit while maintaining relatively low errors on the 
measured fiux. The best fits were obtained by constraining the widths of the doublet 
components to be equal while permitting the relative amplitude of the doublet components 
and the velocity width of the cloud itself to vary. 

The impact of fitting single rows, which subsample the data, was assessed by summing 
3 rows in the spatial direction and then performing the fit. These fits were consistent to 
within 2% for the sum of the [Oil] components. 

We were unable to correct these measurements for extinction because the Ha + [Nil] 
emission line clouds are too confused to extract unblended Ha fiux measurements 
for comparison with the H/3 fiuxes. Our long slit data will reduce the scattered light 
contribution and provide more extended spectral coverage using the Hell lines at 1640 A and 
4686 A for extinction correction. 



6. Results 

Morphology: The NLR emission line clouds of NGC4151 exhibit a biconical distribution. 
This is evident in our STIS [OH] image and the WFPC2 narrow band image (Figure ^ 
and has been illustrated previously with high spatial resolution HST images ( [Evans et 



|al. 19931 , [Boksenberg et al. 1995| , |Winge et al. 1997| , [Hutchings et al. 1998|) . In both the 



dispersed and undispersed images, we find that clouds to the SW of the nucleus are brighter, 
more numerous, and have a more extended distribution than the clouds to the NE. These 
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results are consistent with both ground based (e.g. [Penston et al. 1990| ) and previous HST 



measurements. For [OIII] the emission hne gas in the shtless spectrum is detected ~ 1.3 kpc 
SW of the nucleus and ~ 500 pc to the NE. The [Oil] , H/3, and Ha + [Nil] emission is less 
extended than this, spanning ~ 700 pc SW of the nucleus and ~ 300 pc to the NE (Figures |l] 
and ^). In addition, we detect [01] and [SII] emission line clouds in the NLR with a similarly 
asymmetric spatial distribution. Background subtraction of the STIS dispersed [Oil] image 
(Figure 0) also reveals NLR emission structure of [FeVII] A3760, [Nelll] A3869, Hel A3889, 
and [Nelll] A3967, as well as emission from higher level transitions in the Balmer series (He 
A3970 through H12 A3751). 

Velocities and Dispersions: Our velocity and velocity dispersion measurements of 57 
emission line clouds (Tables 2 and 3; Figures H, ^, and |1^) suggest that there are at least 
two distinct kinematic components present in the NLR of NGC4151. One component 
(LVLVD) is characterized by low velocities and low velocity dispersions (|v| <400kms~^ 
and Av< 130kms~^). A second cloud population (HVHVD) is characterized by high 
velocities and velocity dispersions (|v| >400kms^^ and Av> 130kms^^). An apparent 
intermediate cloud population with low velocities and high velocity dispersions, the LVHVD 
clouds, may be HVHVD clouds whose low velocities are the results of projection effects. 

The HVD clouds all reside within 3.2" (~200pc) of the nucleus. The HVHVD clouds, 
which reside within 1.1" (~70pc), are more tightly clustered about the nucleus than the 
LVHVD clouds. The LVLVD cloud population extends from the nucleus out to a distance 
of ~ 19" (1.2 kpc), with a gap in the cloud distribution at a distance of ~4.0". We may also 
subdivide the LVLVD population into clouds clustered near the nucleus and those beyond 
~4.0" whose dynamics are probably unaffected by nuclear processes. 

Jet-cloud interactions: The high dispersion and highest velocity clouds tend to be 
located near the radio jet axis (Figure |1T]). However, none of the HVHVD clouds are 
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coincident with the radio jet substructure. 



The NGC4151 radio jet contains five distinct knots (Figure |T2|, nomenclature follows 



Uarral et al. 199C) which are not co-aligned along a single position angle, but can be 



represented by an average position angle of 77° ( [Pedlar et al. 19931) . Higher resolution 



measurements indicate that radio knot C4 consists of two components separated by 
92 milliarcseconds (VLBI: [Harrison et al. 1986| and VLBA: [Ulvestad et al. 1998| ). The 



brightest knot within the eastern component is co-spatial with the optical nucleus ([Ulvestad 



et al. 1998| ). Since the radio data presented in Figure |T2| ([Mundell et al. 1995| ) has a 



resolution of 0.16 x 0.15", we aligned the center of the C4 radio component with the optical 
nucleus. To align the two nuclei it was necessary to translate the optical image by 0.06" in 
right ascension and 0.46" in declination. Radio knot C3, at ~0.5" SW of the nucleus, is 
only 0.19" from the projected position of cloud 20. These two inner knots, C4 and C3, lie 
along a position angle of 83°. Outer radio knots C2 and C5, at ~ 0.9" SW and NE of the 



nucleus, lie along a position angle of ~ 73° (Figure |TI]). CI, at ~ 1.9" SW of the nucleus, is 



relatively diffuse and isolated from the [OIII] emission line clouds. 

The NLR [OIII] gas within ~1.5" of the nucleus, a point defined by the 5 a radio 
contour associated with radio knots C2 and C5, is approximately co-spatial with the radio 
emission along PA = 77°. We find that all the clouds in this region are contained within 
an opening angle of ~45° about the NLR axis at 63° with the HVD clouds relatively 
symmetrically distributed about this axis. We have overlaid 5 a radio contours on the 



optical image (Figure [T^ ). Of the 13 clouds contained within the 5cr (1.7mJy) radio 
contour, 4 are LVD clouds. Relaxing the radio contour to 3 a (0.96 mJy), only increases 
the cloud population by 2 HVD clouds. So although the [OIII] emission is approximately 
aligned with the 3 a and 5 o radio contour, the clouds interior to the these contours are 
both LVD and HVD clouds. Within a 1.5" radius of the nucleus but exterior to the 3 a 
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radio contour, we find a similar distribution of LVD to HVD clouds. Of the 10 clouds in 
this region, only 20% are LVD clouds. Thus the radio jet does not appear to be responsible 
for the HVD clouds. In addition, none of the HVHVD clouds are coincident with the radio 
axis and only 5 of 9 HVHVD clouds lie within the 3 a radio contour. 

Cloud 20, whose projected position lies near radio knot C3, is a LVLVD cloud which 
exhibits no kinematic or morphological signatures of interaction with the radio jet. 

Flux Distribution: Flux ratio measurements from different lines in individual clouds 
are consistent with photoionization. 



The radial distribution of H/3 and [OIII] fluxes (Figure |T^ indicates that the flux upper 
envelope decreases as a function of radius for r < 3.2" (~ 200 pc). At larger radii the flux is 
approximately constant. 

Within a ~200pc radius of the nucleus, the [0III]/H/3 flux ratio declines linearly for 
both the redshifted and blueshifted clouds. This behavior, albeit more pronounced for the 
HVD clouds, is exhibited for both the HVD and LVD clouds (Figure p!3|) . This decrease in 
the flux ratio indicates a decreasing ionization parameter ([Ferland and Netzer 1983 ). At 



larger radii the flux ratio is comparable to its value very close to the nucleus. 

To further illustrate the spatial structure in the [0III]/H/3 ratio we have generated a 
[0III]/H/3 image (Figure 113). Much of the discrete cloud structure evident in the dispersed 



[OIII] (Figure |^) and H/3 images is preserved in the ratio image. 

Due to the small number of clouds for which we have [OH] fluxes, we were not able 
to determine if the [0II]/H/3 ratio exhibits the same radial gradient in the inner NLR as 
presented by the high velocity dispersion clouds for [0III]/H/3. 

Continuum NLR Emission: As a further probe of the impact of the radio emission 
on the observed optical emission structure, we have searched the continuum emission 
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associated with the optical clouds near the radio emitting cores for an elevated blue 
continuum flux. This excess flux at shorter wavelengths is a signature of free-free emission 
and scattered light. The continuum flux redward of 3950 A in the [Oil] spectrum and from 
~ 4800 — 5100 A in the [OIII] spectrum was compared with the flux blueward of 3800A . 
Unlike Winge et. al., (1998), we do not see an increased blue continuum for the near nuclear 
clouds in the background subtracted image. Given the complex near-nuclear environment in 
our slitless spectrum and the weak increase in the blue continuum seen in the background 
un-subtracted spectrum it is possible that the background subtraction has masked a weak 
increase in the blue continuum. 

7. Discussion 

As noted in Section |, our measurements (Figure ||) suggest that there are at least two 
distinct kinematic components. One component, the LVLVD clouds, is characterized by 
low velocities (|v| <400kms~^) and low velocity dispersions (Af < 130 km s^^). A second 
cloud population, the HVHVD clouds, appears to be characterized by both high velocities 
and high velocity dispersions (400 < |v| ^1700kms^^, Av >130kms~^). The HVHVD 
cloud population is spatially confined to the inner 1.1" (~70pc) of the NLR, whereas the 
LVLVD clouds are distributed from 0.2" to ~ 19". 

In this section we compare dynamical models with the observed kinematic distribution 
by evaluating the projected radial velocity distribution in the NLR in terms of rotation 
about the nucleus, jet-cloud interactions, and wind driven outflows. Then we discuss the 
relevance of the flux measurements to both dynamical and photoionization models. 

The radial velocity distribution of the NLR clouds (Figure ^ is reminiscent of the 



"S-shape" signature of gravitational rotation about a massive central source (e.g.. Bower 
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et al. 1998). To investigate the significance of tliis distribution, we liave compared tlie 
measured velocity distribution witli tlie velocity profile expected from a simple disk 
undergoing Keplerian rotation. 

The inclination of the model disk was determined by the ratio of the "minor" to 
"major" axis in the inner NLR . This yields an inclination of 70° (Figure ^). In generating 
this model, we have sampled the model at the same positions as the [OIII] velocity 
measurements. Since the [OIII] peak is known to be blueshifted by ~40 kms~^ relative 
to neutral hydrogen, we have displaced the observed velocities by ~ — 40kms^^. Imposing 
this shift improves the agreement between the model and the data (Figure ^. A central 
mass of ~ 4 X lO^M© is required by the model. 

However, as illustrated in Figure |], the data are not well fit by the model. Obviously 
there is a population of high velocity clouds within a 1.1" (~ 70 pc) radius which do not fit 
the model at all. Including these clouds in the model implies a prohibitively high central 
mass of 10^° Mq. In addition, the transition between redshifted and blueshifted velocities 
near the nucleus is too soft to be indicative of gravitational rotation under the influence 
of a central black hole. The model exhibits a sharp transition at the nucleus which is not 
mirrored by the data. Furthermore, the geometry implied by the model requires that the 
axis of rotation of the disk is perpendicular to the radio jet axis, a geometry that would 
be difficult to explain. A recent analysis ([Winge et al. 1999|) of [OIII] A5007 FOG spectra 
augmented with both STIS slitless spectroscopy results ( |Hutchings et al. 1998| ) and ground 



based data interprets the NLR gas kinematics within the central 4" radius as being governed 
by Keplerian rotation about an extended central mass of ~ 1O^M0. However, our analysis 
indicates that the high velocity kinematic component co-existing in this region is not the 
result of simple gravitational motion. 

Other gravitationally dominated models involve orbits which may manifest themselves 
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as bars or triaxial structure. The degree of symmetry across the nucleus requires any of 
these to be rather specially oriented with respect to our line of sight, in order to produce 
the observed radial velocities. In addition, all gravitational models have to be oriented 
such that the observed cloud morphology and velocity structure are not connected to the 
radio axis which is generally considered to be aligned by the angular momentum vector 
of the innermost material around the central mass. We therefore conclude that simple 
gravitational models cannot explain the entire observed velocity distribution within ~ 3" of 
the nucleus. 

At distances exceeding 4" (~ 260 pc) from the nucleus, we find that the gas cloud 
kinematics are consistent with circular rotation in a gravitational field. This is in agreement 
with previous HST and ground based results (e.g. [Fenston et al. 1990|, [Robinson et al. 1994, 



Vila-Vilaro et al. 1995| , [Hutchings et al. 1998| , |Winge et al.T999D . 



As noted in the previous section, we find no detailed correspondence between the optical 
clouds and the radio structure. However, the general distribution of HVD clouds about 
the radio axis may indicate that the radio jets are connected with the driving mechanism 
for these clouds and may be influential in producing their high velocity dispersions. This 
conclusion is consistent with previous results ([Winge et al. 1999| ) which find that the cloud 
kinematics within 4" of the nucleus are disturbed but not dominated by interactions with 
the radio jet. 

The [OIII] flux does not appear to be enhanced for clouds near the radio axis, nor does 
[OIII]/H/? appear to be a function of radio axis location. From our data we see no evidence 
that the jet has affected the ionization state of the NLR gas. 

Since neither gravitational rotation nor jet-cloud interactions are successful in 
describing the observed kinematic distribution of the clouds within ~ 3" of the nucleus, 
another physical mechanism, such as radiation or wind driven outflow, must be responsible 
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for the gas kinematics. In Hutchings et al., (1998), we argued against infall based upon 
opacity and similarity arguments. These arguments still hold. 

In ground based studies with lower spatial resolution, blue asymmetries have been 
observed in the line profiles (e.g. [Heckman et al. 1981| , [Whittle 1985| , Pahari and DeRobertis 



1988|) . Dust extinction within the NLR has been invoked to explain this asymmetry. 
From our observations, the blueshifted clouds SW of the nucleus appear brighter and 
more extended than the redshifted clouds NE of the nucleus. Assuming the SW cone is 
on the far side of the nucleus, our observations are consistent with infalling, intrinsically 
dusty clouds. If the NE cone is on the far side of the nucleus, then our observations are 
consistent with outflowing clouds in a dusty intercloud medium. The outflow scenario 
is compelling for two reasons. First, all the absorbers in the absorption line spectra of 



CIV have outflow velocities ( [Weymann 1997| , Putchings et al. 1998|) . Second, for the high 



velocities in the inner NLR, gravitationally driven infall requires an improbably high central 
mass of ~ 10^° Mq. Furthermore, radio observations ( Pedlar et al. 1993| , [Mundell et al. 



|1995| , pivestad et al. 199^ ) map an outflowing radio jet that is roughly aligned with the 
optical axis of the NLR bicone. 

Emission line ratios of the individual clouds are consistent with photoionization 
models. The [0III]/H/3 ratio declines roughly linearly within an inner ~ 3" (200 pc) radius 
from the nucleus, particularly among the HVD clouds. This decrease in the flux ratio 
indicates a decreasing ionization parameter ( [berland and Netzer 19831 ). Since the ionization 



parameter is proportional to r^^n~^, it appears that the density, n, must decrease as ~r~^. 
At distances further from the nucleus where the clouds have low velocity dispersions, 
the [OIII]/H/? ratio is roughly constant. Within the scatter the [0III]/II/5 ratio in this 
region is comparable to its initial value, indicating that either the ionization parameter is 
roughly constant at larger radii with the density decreasing as r~^ or the data are subject 
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to projection effects, local gradients, or both. 

The ratio of [OI]/[SII] should also be a density indicator. However, the small sample of 
clouds for which we have flux measurements and the uncertainties in these measurements 
could mask a change in the flux ratio as a function of radial distance. The ratios listed 
in Table 5 are consistent with a constant of 0.35 for the 120 pc radius spanned by these 
measurements. 

8. Conclusions 

We have spatially mapped the full velocity field of the NLR for NGC4151. Our data 
support a biconical distribution of emission line gas where clouds to the SW of the nucleus 
are approaching and clouds to the NE have recessional velocities. Within an inner 3.2" 
radius we observe a high velocity dispersion (Av> 130kms~^) population of clouds. The 
[OIII]/H/? flux ratio declines roughly linearly for these clouds. Closer to the nucleus, and 
distributed within a cone of opening angle ~ 45° about the NLR axis, is a population of 
faint HVHVD clouds. 

Some dynamical scenarios have been considered to account for the observationally 
derived kinematic distribution, velocity profiles, and the radial dependence of the flux ratios 
arising in the NLR. Kinematic mapping of the NGC4151 NLR indicates that more than 
one process must be employed to model the entire velocity field associated with the NLR 
emission. Our data are consistent with photoionization by a central source and outflow 
from a spatially unresolved nuclear region. Within an inner 3.2" radius we do not find 
evidence to support gravitationally dominated models such as rotation about a central 
massive source or transport along a nuclear bar structure. 

The evidence we find for radio jet-cloud interactions is provocative but weak. While 
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the HVD clouds are relatively syminetrically distributed within a 45° opening angle cone 
centered on the NLR axis at ~63°, the LVLVD and HVHVD clouds are co-spatial within 
the cones and there is no preference for the HVD or HVHVD clouds to be associated with 
the radio knots which lie along a PA of 77°. 

For the HVD clouds the velocity distribution is sharply peaked about the nucleus and 
the [0III]/H/3 flux ratio declines approximately linearly. 

Our results are consistent with a wind driven outflow as the acceleration mechanism 
for these clouds. Qualitatively the HVHVD clouds resemble the tenuous cloud component 
that is being radially driven outward in NGC 1068 ([Cecil et al. 1990]) . The LVHVD cloud 



population resident within the same projected radius may be participating in the same 
wind driven outflow but its projected velocity along our line of sight is small. Alternate 
driving mechanisms for these clouds are less plausible. Radiative acceleration is inefficient, 
and gravitationally driven infall requires an improbably high central mass of ~ 4 x 10^ Mq 
for radial velocities of 800kms~^. 

Recent HST STIS observations by this collaboration have mapped the NGC 4151 
emission line spectra from 1150 A through 1.0/im at two slit positions which intersect 
multiple clouds in the NLR. We will use these data as input to photoionization models to 
probe the density and ionization structure within the NLR. 

We thank Carole Mundell for providing us with an electronic version of the radio data 
( pVIundell et al. 1995| ) which is presented in Figure |12|. Support for this project was provided 



in response to NASA Announcement of Opportunity OSSA-4-84 through the Hubble Space 
Telescope Project at Goddard Space Flight Center. 
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Fig. 1. — Dispersed G430M image from STIS of the blue region of the spectrum of NGC 4151. 
The positions of the main features in the nuclear spectrum are labelled. The upper panel 
is background subtracted, to enable the faint inner emission clouds to be seen. The lower 
panels show the spectrum as observed, with two stretches to show the inner and outer clouds. 
The wavelength scale is in the observed frame. 

Fig. 2. — Dispersed G430M image from STIS in the observed frame of the H/3 and [O III] 
spectral region of NGC 4151. This may be compared with the image of the same region with 
different spectral dispersion direction, shown in Hutchings et al., 1998. 

Fig. 3. — Dispersed G750M image from STIS of the red region of the spectrum. Two 
stretches are shown to cover the full dynamic range of cloud fluxes. The wavelength scale is 
in the observed frame. 

Fig. 4. — Direct WF image in the light of [O III] (F502N), with the orientation of the later 
observations, to identify the outer clouds measured. 

Fig. 5. — Direct image of inner clouds, labelling them for identification. The image is in the 
light of [O III] taken with the PC and filter F502N. The image is rotated to the orientation 
of the initial STIS dispersed images. The dashed ellipse illustrates the inclination angle used 
in modelling the cloud velocity structure. 

Fig. 6. — Dispersed [O III] 5007 A line image from the STIS, labelling the inner clouds 
measured. The dispersion direction is L-R, and the nuclear and galaxy continuum has been 
subtracted to make individual cloud features more visible. Cloud velocities are given in 
Tables 2 and 3. 

Fig. 7. — This diagram illustrates how cloud locations and velocities are obtained with 
two different roll angles. The thick lines represent the nuclear continuum, and show the 



-29- 



two dispersion directions that occur because of the different roll angles. The images are 
aligned at the nucleus with the known rotational angle differences. In the two dispersed 
images, we see the continuum and emission line cloud images positioned along the same 
dispersion direction, determined by both their spatial location and their radial velocity. The 
undispersed cloud position is located where the lines meet when extended along the two 
dispersion directions. For example, cloud A appears at Al in one image and A2 in the other. 
The location of A is where these two lines intersect, and this is shifted by the same distance 
(in the geometrically corrected image) A-Al and A-A2. The four clouds illustrated show the 
different relative positions they will have with the two roll angles shown. 

Fig. 8. — Radial velocity of clouds with distance from the nucleus perpendicular to the 
dispersion direction (vertical in Figure |: approximately along the radio axis). In the upper 
panel, the two symbols distinguish the broad and narrow lines. The plot is not substantially 
changed if the projected distance from the nucleus is used. The lower panel shows the 
velocity curve for the inclined disk model, with the radial velocities shifted by -40 kms~^, 
as discussed in the text. In the lower panel a model for a rotating thin disk about a central 
mass of 4 X lO^M© is denoted by a dashed line. 

Fig. 9. — Cloud velocity dispersion with distance from the nucleus perpendicular to the 
dispersion direction (vertical distance in Figure^). The high velocity clouds are those which 
have radial velocities in excess of 400 kms~^. 

Fig. 10. — Histogram of cloud velocity dispersions, showing how the bin values change for 
inner subsets of clouds. In all cases the distribution appears bimodal, which we discuss as 
two populations of clouds with dispersions above and below 130 kms~^. 

Fig. 11. — The upper panel illustrates the spatial distribution of the [OIII] emission line 
clouds within 3.2" of the nucleus. The population of HVD, HVHVD, and LVD clouds are 
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denoted, Position angles of the radio cores relative to the nucleus ( [Pedlar et al. 1993| , [Mundell 



et al. 1995| ), the stellar bar ([Pedlar et al. 199^ ), and the extended NLR are illustrated. The 



radio knots C5 and C2, which define a position angle of 73°, are depicted as open squares. 



Fig. 12.— The MERLIN radio map ( [Mundell et al. 1995[ ) overlaid on the F502N PC 



image. While the NLR [OIII] emission is globally aligned with the radio emission, 
there is a lack of detailed correspondence between the [OIII] clouds and the radio knots. 
The resolution of the radio map is 0.15 x 0.16 arcsec. The radio contour levels are at 
[1,2,3,4,5,6,7,8,9,10,15,20,25,30] x 1.7 mJy and 0.96 mJy, where the 5 a contour is at 1.7 mJy 
and the 3cr contour is at 0.96 mJy. 

Fig. 13. — Line fiux measures and ratios for clouds. The symbols distinguish approaching and 
receding clouds in the upper two and lower left panel. In the lower right panel the symbols 
distinguish between high and low velocity dispersion clouds. Note the decline in [0III]/II/9 
in the inner region, which is seen primarily in the lines with high velocity dispersion (lower 
right panel). 

Fig. 14. — Line ratio dispersed image showing [0III]/H/9. Galaxy and scattered nuclear 
continuum light have been removed before deriving the ratio. Dark indicates high ratio, as 
seen in the outer and the high velocity clouds. 
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Table 1. NGC4151 data log 



Date 


Instrument 


Filter 


Grating 


Ac 

[A] 


AA 

[A] 


Exposure 
Time [s 


1995 January 23^ 


WFPC2 


F502N 




5012 


26.8 


860 


1997 March 29^ 


STIS 


Clear 


G430M 


4961 


282 


2690 


1997 July 15 


STIS 


Clear 


G430M 


3843 


284 


3390 


1997 July 15 


STIS 


Clear 


G430M 


4961 


282 


2139 


1997 July 15 


STIS 


Clear 


G750M 


6581 


569 


1960 


1997 July 15 


STIS 


[Oil] 




3740 


80 


3070 


1997 July 15 


WFPC2 


F502N 




5012 


26.8 


1800 


1998 April 5 


WFPC2 


F390N 




3889 


45.3 


1200 


1998 April 5 


WFPC2 


F656N 




6562 


22 


700 


1998 April 5 


WFPC2 


F658N 




6590 


28.5 


402 


1998 April 5 


WFPC2 


F673N 




6733 


47.2 


1200 



^Archival data 
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Table 2. Cloud Velocities 



Cloud # Radial Velocity^ a^v At; (Dispersion) ctav 

kms~^ kms~^ kms^^ kms^^ 



1 
2 
3 

4 

5 

6 

7 

8 

9 

9b 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

19a 

20 

21^ 



230 


24 


47 


25 


66 


40 


278 


14 


192 


12 


11 


13 


105 


- 


265 


18 


72 


10 


10 


16 


495 


21 


36 


13 


337 


36 


190 


14 


125 


21 


-93 


1 


•152 




•107 


17 


235 


47 


•122 


29 


•215 


- 


36 


19 


•538 


10 



183 


74 


340 


14 


63 


13 


44 


27 


75 


5 


153 


38 


94 


- 


181 


24 


83 


20 


241 


4 


160 


45 


223 


7 


300 


14 


24 


21 


205 


6 


282 


36 


317 


115 


248 


128 


157 


38 


396 


43 


79 


- 


78 


12 


294 


_ 
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Table 2 — Continued 



Cloud # Radial Velocity^ ajiy A t> (Dispersion) a/^^ 

kms~^ kms"^ kms~^ knis~^ 



21x'= 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 



-43 


- 


-279 


19 


-280 


24 


NC^ 




-248 


16 


-814 


53 


-161 


29 


190 


32 


84 


21 


78 


15 


20 


- 


;host 




7 


21 


7 


26 


7 


14 


-3 


8 


24 


13 


27 


14 


6 


25 





1 


8 


18 


-62 


9 


-7 


_ 



116 


6 


198 


45 


243 


10 


304 


4 


119 


169 


156 


141 


70 


71 


13 


14 


210 


- 


64 


55 


43 


16 


48 


- 


76 


25 


34 


9 


46 


10 


41 


37 


51 


- 


59 


51 


68 


63 


74 


_ 
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Table 2 — Continued 



Cloud # Radial Velocity^ aRv A v (Dispersion) cta,; 

kms^^ kms^^ kms^^ kms^^ 

44 5 - - 

45 -18 - - 

46 -30 - - 

47 -5 - 45 - 

48 4 - - 

49 -13 - 101 

50 -20 - - 

51 -20 - - 



^Relative to the heliocentric recessional velocity of 997 km s ^ ( [Pedlar 
et al. 19931 ) 



b 



Velocity is from earlier epoch ( [Hutchings et al. 19981 ), this cloud is 



blended with nucleus in this epoch. 

'^Closest UNblended cloud to the nucleus 

'^Not Confirmed 
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Table 3. High Velocity Clouds 



Cloud # RV Uncertainty Ry a^^ dxo 
kms"^ kms~^ kms~^ 



dyo Neighboring Cloud 



a 


-1716 


256 


175 


0.493 


-0.963 


19 


b 


-1516 


261 


140 


0.028 


-0.963 


22 


c 


-1107 


187 


278 


0.231 


-0.913 


19, 19a 


d 


-1448 


248 


291 


0.606 


-0.659 


15 


n 


-1162 


412 


393 


0.078 


-0.254 


21x, 20 


g 


846 


104 


265 


-0.253 


0.811 


8 
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Table 4. H/9 and [OIII] Cloud Fluxes Derived From Gaussian Fits 



Cloud # Rp Flux [OIII] Flux [O III]/H/3 

erg s^^ cm^^ erg s^^ cm~^ 



1 3.09x10"^^ 3.29x10-13 xo.65 

2 1.35x10-1^ 1.42x10-13 ^0.52 

3 8.91x10-16 9.71x10-1^ 10.90 

4 4.92x10-15 4.90x10-1^ 9.96 

5 2.76x10-15 3.14x10-1^ 11.38 
6+11 2.92x10-1^ 3.37x10-13 11.54 

7 3.45x10-15 3.45x10-1^ 10.00 

8 8.33x10-15 9.77x10-1^ 11.73 
9a 1.93x10-1^ 2.13x10-13 11.04 
9b 2.30x10-1^ 2.45x10-13 10.65 
10 1.77x10-1^ 1.57x10-13 8.87 

12 1.30x10-1^ 1.48x10-13 11.38 

13 2.38x10-15 2.12x10-1^ 8.93 

14 3.02x10-15 2.17x10-1^ 7.19 
15+17 2.77x10-1^ 3.41x10-13 12.31 

16 1.09x10-1^ 1.25x10-13 11.47 

18 1.06x10-1^ 1.02x10-13 9.62 

19 4.58x10-1^ 4.60x10-13 10.04 
19a 1.43x10-1^ 1.42x10-13 9.93 

20 4.67x10-1^ 5.84x10-13 12.51 
21x 2.64x10-1^ 2.91x10-13 11.02 

22 2.28x10-1^ 3.07x10-13 13.46 

23 3.27x10-1^ 3.53x10-13 10.80 
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Table 4 — Continued 



Cloud # H/5 Flux [OIII] Flux [O III]/H/3 





erg s ^ cm ^ 


erg s 1 cm ^ 




25 


1.72x10"^^ 


1.46x10-13 


8.49 


26 


7.11x10-15 


1.21x10-13 


17.02 


27 


6.17x10-15 


5.78x10-1^ 


9.37 


28 


2.90x10-15 


2.03x10-1^ 


7.00 


29 


7.89x10-16 


1.19x10-1^ 


15.08 


30 


9.43x10-1"^ 


1.31x10-1^ 


13.89 


35 


3.98x10-16 


5.70x10-15 


14.32 


37 


5.90x10-16 


8.26x10-15 


14.00 


39 


3.89x10-16 


4.54x10-15 


11.67 


41 


2.93x10-16 


3.12x10-15 


10.65 


42 


5.20x10-16 


6.86x10-15 


13.19 
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Table 5. NGC4151 Cloud Fluxes Derived From Gaussian Fits (xlO ^^ erg s ^ cm 



-2 



Cloud 




[Oil] 


H/3 


[OIII] 


[OI] 


[SII] 


[OII]/H,3 


[Oni]/H/3 




Flux 


Error 


Flux 


Error 


Flux 


Error 


Flux 


Error 


Flux 


Error 






1 


4.11 


0.49 


2.60 


0.69 


26.9 


6.32 


1.36 


0.58 


2.51 


0.21 


1.58 


10.35 


6+11 


2.93 


0.91 


1.55 


0.33 


19.3 


4.22 


1.18 


0.16 


3.08 


0.51 


1.89 


12.45 


9a 


5.83 


1.78 


1.93 


0.54 


21.3 


6.17 


1.23 


0.33 


3.41 


0.50 


3.02 


11.04 


15 


5.77 


1.33 


1.18 


0.42 


13.2 


4.42 


2.89 


0.52 


5.35 


-0.06 


4.89 


11.19 


19 


6.43 


2.28 


3.77 


0.85 


41.20 


6.36 


2.89 


0.99 


5.00 


2.19 


1.71 


10.93 


20 


10.00 


2.90 


4.09 


1.26 


51.70 


5.17 


3.33 


1.06 


9.96 


2.24 


2.44 


12.64 


22 


8.59 


1.65 


1.88 


0.63 


23.50 


9.40 


0.83 


0.09 


4.74 


-0.41 


4.57 


12.50 


23 


6.92 


1.99 


3.27 


1.04 


35.30 


8.95 


2.10 


0.47 


5.82 


0.97 


2.12 


10.80 


25 


4.14 


1.24 


1.24 


0.45 


10.80 


3.96 


0.61 


0.02 


2.73 


0.76 


3.34 


8.71 


30 


0.55 


0.12 


0.12 


0.04 


1.65 


0.53 










4.62 


13.87 
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